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SUMMARY 

The  aim  of  this  report  is  to  present  a  method  for  predicting 
the  aerodynamic  characteristics  of  slender  bodies  of  revolution 
at  small  incidence,  under  flow  conditions  such  that  the  boundary 
layer  is  turbulent. 

Firstly  a  panel  method  based  on  slender  body  theory  is 
developed  and  used  to  calculate  the  surface  velocity  distribution 
on  the  body  at  zero  incidence.  Secondly  this  velocity 
distribution  is  used  in  conjunction  with  an  existing  boundary 
layer  estimation  method  to  calculate  the  growth  of  boundary 
layer  displacement  thickness  which  is  added  to  the  body  to 
produce  the  effective  aerodynamic  profile.  Finally,  recourse 
is  again  made  to  slender  body  theory  to  calculate  the  normal 
force  curve  slope  and  centre  of  pressure  position  of  the 
effective  aerodynamic  profile. 

Comparisons  made  between  predictions  and  experiment  for  a 
number  of  slender  bodies  extending  from  highly  boattailed 
configurations  to  ogive-cylinders,  and  covering  a  large  range 
of  boundary  layer  growth  rates,  indicate  that  the  method  is 
useful  for  missile  design  purposes. 
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Secondly. this  velocity  distribution  is  used  in  conjunction  with  an  existing 
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displacement  thickness  which  is  added  to  the  body  to  produce  the  effective 
aerodynamic  profile.  Finally,  recourse  is  again  made  to  slender  body  theory 
to  calculate  the  normal  force  curve  slope  and  centre  of  pressure  position  of  the 
effective  aerodynamic  profile. 

Comparisons  made  between  predictions  and  experiment  for  a  number  of  slender 
bodies  extending  from  highly  boattailed  configurations  to  ogive-cylinders,  and 
covering  a  large  range  of  boundary  layer  growth  rates,  indicate  that  the  method 
is  useful  for  missile  design  purposes. 
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1.  INTRODUCTION 


WSRL-0093-TR 


Calculations  based  on  potential  flow  aerodynamic  theories  generally  provide 
reasonably  close  approximations  to  the  lift  and  pitching  moment  experienced  by 
near-planar  configurations  such  as  aerofoils  and  wings  at  small  angles  of  inci¬ 
dence.  However,  similar  calculations  applied  to  extreme  three-dimensional 
shapes  such  as  bodies  of  revolution  often  lead  to  force  and  moment  estimates 
which  are  grossly  in  error,  particularly  if  the  body  has  a  boattail.  The  main 
source  of  error  is  that  potential  flow  methods  do  not  admit  the  development  of  a 
boundary  layer,  which,  at  the  Reynolds  numbers  of  practical  interest,  is  predom¬ 
inantly  turbulent.  For  an  aerofoil  or  wing  the  boundary  layer  displacement 
thickness  provides  a  thin  sheath  about  the  wing  surface,  and  this  causes  a  small 
perturbation  to  lift  and  pitching  moment.  Wing  lift  calculations  are  often 
refined  by  allowing  for  the  displacement  thickness  of  the  boundary  layer  as  an 
additional  contribution  to  the  aerofoil  thickness  distribution.  For  a  body  of 
revolution  at  small  incidence  the  boundary  layer  thickness  can  grow  very  rapidly 
(particularly  if  the  body  has  a  boattail,  in  which  case  the  pressure  gradient  is 
adverse  over  the  rear  of  the  body,  and  at  the  same  time  the  streamlines  are 
converging)  and  the  displacement  thickness  of  the  boundary  layer  envelops  the 
body  in  a  sheath  which  produces  a  substantially  changed  aerodynamic  profile,  with 
consequent  large  changes  in  aerodynamic  properties. 

This  feature  can  be  readily  illustrated  using  the  results  of  slender  body 
theory.  Thus  (see,  for  example.  Ward  (1))  the  normal  force  (-Z)  and  pitching 
moment  (M)  about  the  nose  on  a  slender  body  of  revolution  with  pointed  nose  and 
incidence  u,  are  given  by 


-Z  =  PUo2a  S(L) 


M  =  -pU02  a  {  LS(L)  -V}  J 

where  L  is  the  body  length, 

S(L)  is  the  base  cross-sectional  area, 

Uo  is  the  free  stream  speed, 

V  is  the  volume  of  the  body,  and 
p  is  the  fluid  density. 

When  the  body  has  a  boattail,  S(L)  can  be  much  less  than  the  maximum  cross- 
sectional  area,  and  equation  (1)  shows  that  Z  will  then  be  small  and  M  is  likely 
to  be  positive.  If  slender  body  theory  is  applied  to  the  profile  formed  by  the 
addition  of  the  boundary  layer  sheath  (the  displacement  thickness)  to  the  body 
profile,  then  the  effective  base  area  will  be  substantially  greater  than  S(L) 
and  V  will  be  increased.  The  net  effect  is  to  cause  a  substantial  increase  in 
normal  force  (perhaps  by  a  factor  of  2)  and  in  most  cases  a  very  much  less 
positive  pitching  moment.  The  effect  of  the  boundary  layer  displacement  thick¬ 
ness  is  thus  seen  to  be  very  significant,  and  should  be  included  in  aerodynamic 
estimates,  particularly  when  boattailed  bodies  are  being  examined. 

This  is  not  a  new  idea  and  has  been  in  vogue  for  some  decades.  For  example, 
it  was  used  by  Young  (2)  and  Myring  (3)  in  estimating  the  drag  of  bodies  of 
revolution,  while  its  use  is  basic  to  theoretical  estimates  of  the  Magnus  force 
and  moment  experienced  by  slender  rotating  bodies  of  revolution  (see  for  example, 
references  4  to  7) .  Recently,  Esch  (8)  included  the  boundary  layer  displacement 
thickness  in  his  comprehensive  investigation  into  the  effects  of  Reynolds  number 
and  transition  position  on  the  normal  force  and  pitching  moment  characteristics 
of  long  slender  ogive-cylinders. 
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The  aim  of  the  present  study  has  been  to  develop  a  method  of  calculating  the 
normal  force  and  centre  of  pressure  position  of  bodies  of  revolution  at  subsonic 
speeds  and  low  incidence,  covering  all  slender  shapes  irrespective  of  whether 
they  have  boattails  or  not.  Simplicity  has  been  a  prime  consideration  in  the 
development  of  each  step  of  the  method,  even  at  the  expense  of  some  loss  in 
accuracy.  At  the  outset  it  is  stated  that  a  universal,  accurate,  simple  method 
was  not  achieved,  but,  provided  that  the  body  is  reasonably  slender,  estimates 
are  obtained  which  are  adequate  for  many  design  purposes. 

The  method  is  divided  into  three  parts.  First  the  inviscid-flow  velocity 
distribution  on  the  body  surface  is  calculated;  then  the  distribution  of  boun¬ 
dary  layer  displacement  thickness,  which  defines  the  aerodynamic  sheath  surround¬ 
ing  the  body  and  thus  the  effective  aerodynamic  profile;  and  finally  the  normal 
force  and  centre  of  pressure  position  are  calculated  by  applying  slender  body 
theory  to  the  effective  aerodynamic  profile.  In  Section  2  a  panel  method  based 
on  slender  body  theory  is  outlined  (and  developed  in  Appendix  I)  to  enable 
calculation  of  the  velocity  distribution.  The  turbulent  boundary  layer  method 
of  Stratford  and  Beavers  (9)  is  adapted  in  Section  3  for  calculating  the  displace¬ 
ment  thickness,  while  in  Section  4  these  methods  are  applied  to  a  number  of 
configurations,  and  the  estimated  normal  force  characteristics  and  centre  of 
pressure  positions  are  compared  with  experimentally  determined  values. 


2.  VELOCITY  DISTRIBUTION 

The  pressure  gradients  over  the  nose  and  the  downstream  parallel  section  of  a 
body  of  revolution  are  such  that  the  boundary  layer  growth  rate  is  small.  In 
these  regions  the  aerodynamic  profile  is  therefore  little  different  in  shape  from 
the  body  profile,  and  errors  in  calculated  boundary  layer  thickness  distribution 
arising  from  errors  in  surface  velocity  distribution  are  unimportant.  Hence 
over  the  nose  and  parallel  section  some  simplifying  approximations  can  be  used 
in  calculating  the  velocity  distribution.  In  the  vicinity  of  the  boattail  the 
boundary  layer  grows  very  rapidly,  due  to  the  combined  effect  of  an  adverse 
pressure  gradient  and  also  to  the  convergence  of  streamlines,  which  provides  a 
mechanism  by  which  the  thickening  boundary  layer  is  concentrated  and  thus  thickens 
further.  It  will  be  demonstrated  in  Section  4  that  even  here  some  errors  in 
velocity  distribution  can  be  tolerated.  It  appears  that  the  rate  of  boundary 
layer  growth  caused  by  streamline  convergence  is  so  significant  that,  in  compari¬ 
son,  errors  introduced  by  moderate  inaccuracies  in  the  representation  of  velocity 
gradient  are  relatively  unimportant. 

Hence,  summing  up,  the  boundary  layer  behaviour  over  all  sections  of  the  body 
can  be  adequately  predicted  using  an  approximate  representation  to  the  velocity 
distribution. 

As  discussed  in  Section  3  the  boundary  layer  displacement  thickness  distribu¬ 
tion  is  calculated  at  zero  incidence  only  and  is  added  to  the  body  to  find  an 
effective  aerodynamic  profile  which  is  taken  to  hold  for  all  (small)  incidences. 

It  is  therefore  only  necessary  to  calculate  the  velocity  distribution  at  zero 
incidence,  and  the  method  selected  is  based  on  an  application  of  slender  body 
theory  by  Spreiter  and  Alksne  (10),  who  derived  an  expression  particularly 
suitable  for  this  purpose.  Thus,  considering  the  uniform  flow  at  speed  Uo  over 
a  smooth  pointed  slender  body  of  revolution  at  zero  incidence,  the  perturbation 
velocity  u  is  given  by 


u^ 

Uo 


S"M 

2jt 


In 


fir  _ 

V  x(L-x) 


(2) 


where  L  is  the  body  length, 

S  is  the  body  cross-sectional  area  distribution, 


p  =  /mt  ,  Nfe  being  the  freestream  Mach  number. 


and  £  is  a  dummy  variable. 

Equation  (2)  becomes  particularly  simple  to  evaluate  if  the  S"  terms  are 
either  constant  or  vary  linearly  with  x.  Now,  most  missile  body  configurations 
are  composed  of  a  conical  or  ogival  nose, one  or  more  parallel  sections,  and 
perhaps  a  boattail  or  flare.  S"(x)  is  constant  for  a  conical  nose,  a  boattail 
and  a  flare,  and  is  zero  for  a  parallel  section.  Furthermore,  a  linear  distri¬ 
bution  of  S"(x)  of  the  form  S"(x)  =  ax  +b  may  be  shown  (see  Appendix  I)  to 
represent  a  curved  nose  profile  which,  in  the  context  of  this  investigation,  can 
be  taken  as  an  approximate  representation  of  all  ogival  nose  profiles.  Hence, 
various  sections  of  a  body  of  revolution  can  be  considered  as  panels,  in  each  of 
which  S"(x)  is  either  constant  or  varies  linearly.  If  equation  (2)  is  solved 
for  these  cases  the  total  perturbation  velocity  can  be  found  by  summing  together 
the  contributions  produced  by  the  individual  panels.  It  will  be  appreciated 
that  the  representation  of  a  body  by  a  number  of  panels  in  the  manner  described 
above  is  likely  to  lead  to  a  discontinuous  distribution  of  S"(x).  At  each 
discontinuity  slender  body  theory  breaks  down  and  indicates  a  singular  behaviour 
in  the  perturbation  velocity.  However,  it  usually  happens  that  S"(x)  discontinui¬ 
ties  are  located  sufficiently  far  apart  for  the  regions  of  singular  behaviour  to 
be  reasonably  separated,  and  they  can  be  smoothed  out  and  still  permit  an  adequate 
prediction  of  the  velocity  distribution.  Examples  to  be  given  in  Section  4  will 
illustrate  this  point. 

The  three  S"(x)  distributions  illustrated  in  figure  1  require  examination  and 
the  appropriate  perturbation  velocities  are  evaluated  in  Appendix  1,  which  also 
includes  the  application  of  S"(x)  distributions  to  noses,  boattails  and  stabilis¬ 
ing  flares. 
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3.  BOUNDARY  LAYER  DEVELOPMENT 

When  a  body  of  revolution  is  at  incidence  a  different  boundary  layer  displace¬ 
ment  thickness  distribution  develops  along  each  meridian,  and  the  effective 
aerodynamic  body  will  not  be  a  body  of  revolution.  If,  as  mentioned  in  Section 
1,  this  meridianal  variation  is  ignored  and  the  boundary  layer  growth  is  calcu¬ 
lated  at  zero  incidence  only,  the  method  of  approach  is  kept  simple,  but  at  the 
expense  of  some  loss  in  accuracy.  In  an  endeavour  to  estimate  the  magnitude  of 
this  loss  an  investigation  was  conducted  in  which  for  one  slender  boattailed  body 
(body  no.  1)  the  velocity  distribution  was  calculated  on  the  extreme  windward  and 
leeward  meridians  at  5°  incidence.  The  boundary  layer  displacement  thickness 
distributions  were  then  estimated  and  were  found  to  differ  considerably,  their 
mean  being  close  to  the  distribution  calculated  on  the  body  in  axisymmetric  flow. 
However,  in  this  report  the  effective  aerodynamic  profiles  rather  than  boundary 
layer  thickness  distributions  should  be  compared.  Such  a  comparison  revealed  a 
small  relative  difference  between  the  effective  aerodynamic  profiles  at  5°  and 
at  0°  incidence. 

A  further  effect  of  different  boundary  layer  growth  rates  on  the  windward  and 
leeward  sides  is  that  the  effective  aerodynamic  profile  is  not  symmetrical  but 
has  a  curved  centreline,  which  of  course  influences  the  aerodynamic  loading. 

For  the  particular  case  investigated  above  the  non-uniform  boundary  layer  growth 
yielded  an  effective  aerodynamic  boattail  with  an  incidence  about  3%  less  than 
the  geometrical  value,  which  was  calculated  (using  an  application  of  slender  body 
theory  reported  by  Isaacs  (11))  to  give  a  5%  increase  in  normal  force  curve  slope 
and  a  forward  centre  of  pressure  movement  of  less  than  0.3  calibres.  Although 
these  increments  are  significant  (particularly  that  to  normal  force  curve  slope) 
the  method  proposed  does  not  include  an  allowance  specifically  to  cover  boundary 
layer  induced  camber.  The  boundary  layer  displacement  thickness  distribution 
on  the  body  at  zero  incidence  is  calculated  and  added  to  the  body  to  find  an 
effective  aerodynamic  profile  that  is  taken  to  hold  for  all  (small)  incidences. 

As  will  be  described  shortly  an  empirically  determined  value  of  the  boundary 
layer  shape  factor  is  used  in  the  boattail  region  to  enable  calculation  of  the 
displacement  thickness  there.  The  choice  of  shape  factor  is  such  that  in  most 
cases  the  error  in  estimated  normal  force  curve  slope  is  considerably  less  than 
5%. 

Stratford  and  Beavers  (9)  formulated  a  convenient  method  for  calculating  the 
turbulent  boundary  layer  properties  of  two  and  three-dimensional  shapes  in 
compressible  flow.  At  any  point  the  boundary  layer  is  represented  as  that 
growing  on  an  equivalent  flat  plate  at  a  distance  from  the  leading  edge  dependent 
upon  surface  Mach  number  distribution  and  body  shape.  Now,  slender  body  theory 
shows  that  compressibility  effects  are  not  very  pronounced  for  slender  bodies  of 
revolution  (as  shown  by  the  weak  influence  of  the  terms  in  P  in  equations  (1.1) 
(1.2)  and  (1.3)  and  therefore  an  incompressible  version  of  Stratford  and  Beavers' 
method  is  considered.  Thus,  the  boundary  layer  at  distance  x  from  the  body  nose 
is  the  same  as  that  on  a  flat  plate  at  distance  X  from  the  leading  edge  given  by 


where  U  is  the  fluid  speed  at  the  body  surface, 
7  is  a  constant  such  that 
7  -  1.25  when  Rx  ~  106  , 

7  =  1.20  when  Rx  ~  107  , 

UX 

X  3  v 


and  R. 
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is  a  Reynolds  number  based  on  local  stream  speed  outside  the  boundary  layer,  v 
being  the  kinematic  viscosity. 

The  momentum  thickness  0  is  given  by 


e 

6 


0.036  X  Rx  for  Rx  ~  106  , 


0.022  X  Rx  for  Rx  ~  107  . 


(3) 


Stratford  and  Beavers  recommend  determination  of  the  displacement  thickness  5* 
for  small  pressure  gradients  by  assuming  a  l/7th  power  law  velocity  profile,  from 

6  * 

which  the  shape  factor  H  =  5-  can  be  evaluated  and  thence  6*  using  equation  (3). 
Thus, 

-5 

8*  =  0.036  HXRX  for  Rx  ~  106  , 


8*  =  0.022  HXR^  for  Rx  ~  107  , 


For  severe  pressure  gradients,  Stratford  and  Beavers  recommend  calculating  the 
variation  of  H  over  the  body  surface,  and  thence  8*  from  (4).  The  method  adopted 
in  this  report  is  to  treat  the  body  in  two  sections,  namely  that  part  over  which 
the  cross-sectional  area  (or  radius)  increases  with  axial  distance  or  is  constant, 
and  that  part  over  which  it  decreases  with  axial  distance  (i.e.,  the  boattail 
section).  Over  the  first  part  the  boundary  layer  growth  is  generally  similar  to 
that  over  a  flat  plate,  for  which  the  commonly  accepted  value  of  H  =  1.4  is  used. 
Over  a  boattail  the  pressure  gradient  can  be  very  severe  indeed  and  H  will 
increase  substantially  beyond  1.4.  Reference  12  contains  several  sets  of  experi¬ 
mental  data  relevant  to  the  flow  over  bodies,  and  examples*  roughly  representative 
of  boattail  flows  give  values  of  H  between  2  and  4. 

To  simplify  boundary  layer  calculations  over  the  boattail,  the  variation  of  H 
with  axial  distance  from  1.4  to  perhaps  greater  than  3  has  been  replaced  by  a 
constant  value  for  all  boattails,  chosen  empirically  using  the  following  approach. 
Equation  (1)  shows  that  the  body  normal  force  depends  on  the  base  area  S(L)  of  the 
effective  aerodynamic  profile,  and  that  S(L)  is  also  a  dominant  parameter  in 
determining  the  pitching  moment  and  therefore  the  centre  of  pressure  position. 


Now 


S(L) 


where  subscript  B  refers  to  conditions  at  the  base.  Any  errors  in  8g  will  be 

of  greatest  significance  on  configurations  for  which  rfi  is  small,  i.e.,  on 

body  shapes  with  extensive  boattailing.  It  follows  that  if  H  is  chosen  empiri¬ 
cally  it  should  be  selected  in  such  a  way  that  the  predicted  and  measured  normal 
force  and  centre  of  pressure  positions  match  as  closely  as  possible  for  bodies 
with  extensive  boattails,  since  these  are  the  key  configurations.  If  this 
chosen  H  is  then  applied  to  evaluate  8*  on  a  body  with  less  severe  boattailing, 

it  will  overestimate  8*  ,  but  because  r^  is  now  not  so  small  the  errors  in  S(L) 

*  Relevant  tramples  are  those  numbered  3600,  3700,  3800,  4000,  5000  and  5100  in 
reference  12. 
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(and  therefore  in  normal  force  and  pitching  moment)  are  not  likely  to  be  very 
significant.  Thus,  a  single  H  value  can  be  chosen  to  cover  all  boattails. 

By  examining  several  boattailed  bodies  it  was  found  that  H  =  3.0  led  to  the  best 
overall  estimate  of  normal  force  and  centre  of  pressure  position  for  boattailed 
bodies . 

Summing  up,  over  that  part  of  a  body  where  the  radius  either  increases  with 
axial  distance  or  is  constant,  H  is  taken  to  be  1.4.  When  the  radius  decreases 
with  axial  distance  H  is  taken  to  be  3.0.  The  appropriate  value  of  H  is 
substituted  into  equation  (4)  for  the  evaluation  of  6*. 

This  simple  choice  of  two  constant  values  for  H  can  lead  to  an  uncertainty 
when  considering  a  configuration  with  a  parallel  afterbody.  Thus,  the  after¬ 
body  can  be  thought  of  either  as  a  section  having  constant  radius,  in  which  case 
H  =  1.4,  or  as  the  limiting  case  for  a  boattail,  namely  a  zero-angle  boattail, 
in  which  case  H  =  3.0.  These  two  alternatives  have  been  examined  in  Section  4 
in  connection  with  two  ogive-cylinders  (bodies  no.  5  and  6),  and  the  differences 
were  found  to  be  small.  This  finding  is  consistent  with  the  statement  made  above 


B  * 

that  for  such  bodies  —  is  small  and  changes  due  to  errors  in  estimates  of  6R 

rB 


B 


are  not  likely  to  be  very  significant. 

It  will  be  noted  that  since  H  changes  its  value  discontinuously  at  the  junc¬ 
tion  of  a  parallel  section  and  a  boattail,  5  will  change  discontinuously  there 
and  the  effective  aerodynamic  profile  will  also  have  a  discontinuity.  This 
reflects  the  crudeness  of  the  method  in  replacing  a  varying  shape  factor  over 
the  boattail  by  a  constant  value.  Even  though  slender  body  theory  does  not  hold 
in  the  vicinity  of  this  discontinuity,  nevertheless  equation  (1)  can  be  applied 
to  evaluate  both  normal  force  and  pitching  moment  and  thence  the  centre  of 
pressure  position  of  the  effective  aerodynamic  profile. 


4 .  EXAMPLES 

Several  examples  are  given  below  in  order  to  illustrate  how  well  the  surface 
velocity  distribution  can  be  predicted  by  remarkably  few  panels  of  either 
constant  or  linearly  varying  S"(x),  and  also  how  closely  the  method  predicts  the 
normal  force  and  centre  of  pressure  position.  Although  experimental  normal 
force  and  centre  of  pressure  position  data  are  presented  at  speeds  high  enough 
for  the  effects  of  compressibility  to  be  experienced,  in  the  interests  of 
simplicity  all  velocity  distributions  and  boundary  layer  calculations  have  been 
based  on  incompressible  flow  considerations. 

4.1  Slender  boattailed  bodies 

4.1.1  Body  no .  1 

The  slender  bomb  body  configuration  (figure  2(a))  was  represented 
by  five  panels  of  S"(x)  distribution,  as  shown  in  figure  2(b).  The 
curved  nose  was  represented  by  the  "universal"  nose  profile  (equation 
(1-4))  and  the  corresponding  S"(x)  distributions  comprising  two  panels, 
the  boattail  by  a  single  panel  having  a  constant  S"(x)  distribution, 
and  two  panels  representing  the  curved  transition  region  between  the 
parallel  centre  section  and  the  boattail.  The  crosses  in  figure  2(a) 
show  the  profile  produced  by  these  panels  and  the  agreement  with  the 
actual  profile  is  seen  to  be  close.  The  deduced  velocity  distribution 
is  given  in  figure  2(c).  Also  shown  in  this  figure  is  a  velocity 
distribution  calculated  by  Haselgrove  (13)  who  adapted  a  computer 
programme  of  Albone  (14),  results  of  which  are  highly  accurate  for 
slender  shapes.  Apart  from  singular  regions  occurring  at  the  beginn¬ 
ing  and  end  of  S"(x)  distributions,  the  two  velocity  distributions 
agree  closely.  Furthermore,  except  near  the  nose  point,  the  most 
obvious  line  through  the  singular  regions  provides  a  close  estimate  of 
the  velocity  distribution. 


IhI-J 


Figure  2.  Approximation  to  body  no.  1  and  calculated  velocity 

distribution  at  M  =  0,  a  =  0° 
o  * 

(a)  Body  profile 

(b)  S" (x)  distribution 


(c)  Velocity  distribution 
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The  discrepancy  near  the  nose  point  is  significant  but  it  will  not  lead 
to  large  errors  in  estimates  of  boundary  layer  displacement  thickness 
distribution  because  of  the  highly  favourable  pressure  gradient  and 
the  expanding  body  surface  with  axial  distance,  both  of  which  lead  to 
a  thin  boundary  layer.  The  comparison  in  figure  2(c)  provides  an 
indication  that,  at  zero  incidence  at  least,  the  body  in  figure  2(a) 
behaves  as  if  slender  everywhere  except  near  the  nose. 

The  shape  in  figure  2(a)  was  tested  in  a  wind  tunnel  by  Landers  (15) 
who,  in  order  to  allow  for  an  internal  strain  gauge  force  measurement 
balance,  had  to  distort  the  rear  of  the  body  as  shown  in  the  region  A 
in  figure  2(a).  In  these  experiments  Landers  used  a  nose  boundary 
layer  trip  to  ensure  that  a  turbulent  boundary  layer  occurred;  he  also 
checked  that  the  flow  was  attached  over  the  boattail.  The  experimental 
results  of  normal  force  coefficient  slope  and  centre  of  pressure  posi¬ 
tion  are  given  in  Table  1.  Also  shown  are  the  results  of  slender  body 
theory  assuming  inviscid  flow,  and  in  addition  the  estimated  values 
based  on  the  method  of  this  paper,  using  the  velocity  distribution  in 
figure  2(c).  For  both  of  these  estimates,  allowance  was  made  for  the 
base  distortion  shown  in  figure  2(a).  The  present  method  is  seen  to 
provide  a  large  correction  to  slender  body  theory  and  to  predict  normal 
force  and  centre  of  pressure  characteristics  which  are  in  reasonably 
close  agreement  with  experiment. 

4.1.2  Body  no.  2 

The  body  profile  in  figure  3(a)  was  represented  approximately  by 
four  panels  (figure  3(b)),  two  representing  the  nose  and  two  the 
boattail.  The  profile  thus  produced  is  seen  to  match  the  body  closely 
except  for  the  front  half  of  the  nose  which  is  too  fat. 


(a)  Body  profile 

(b)  S"(x)  distribution 


Figure  3.  Approximation  to  body  no.  2 


Calculated  slender  body  characteristics  and  those  given  by  the  present 
method  are  shown  at  the  second  entry  in  Table  1,  where  they  are  compared 
with  the  measurements  of  Marsden  (16)  made  at  high  Reynolds  number. 

Again  the  defects  of  slender  body  theory  in  inviscid  flow  are  revealed, 
as  well  as  the  gratifying  agreement  produced  by  the  present  method. 

4.2  Non-slender  boattailed  bodies 

The  next  two  examples  concerning  the  two  bodies  shown  in  figure  4  are 
interesting  in  that  they  indicate  what  might  happen  when  the  method  is 
applied  to  configurations  which  do  not  satisfy  the  slenderness  requirements 
of  slender  body  theory.  Both  bodies  have  the  same  blunt  nose  and  a  relative¬ 
ly  steep  boattail  angle.  Body  no.  3  boattail  extends  to  the  base  while  no.  4 

has  a  parallel  afterbody.  A  close  match  to  both  shapes  was  achieved  by  using 
two  panels  to  represent  the  nose  and  one  to  represent  the  boattail  (S"(x) 
distribution  A  in  figure  4(b)).  In  addition,  further  panels  were  also  used 
to  fair  the  sharp  kinks  in  body  profile  at  either  end  of  the  boattail,  thus 
weakening  the  singularities  in  their  vicinity.  The  velocity  distributions 
calculated  from  the  panels  are  shown  in  figure  4(c)  which  also  contains  a 
much  more  accurate  velocity  distribution  calculated  by  Haselgrove  (13) .  In 
spite  of  the  close  match  to  the  body  profiles  illustrated  in  figure  4(a), 
slender  body  theory  produces  substantial  errors  in  velocity  distribution, 
particularly  in  the  vicinity  of  the  nose.  The  addition  of  extra  panels  to 
reduce  singular  behaviour  at  either  end  of  the  boattail  is  seen  to  produce 
significant  improvements  in  velocity  distribution.  Even  so,  the  singular 
regions  are  still  extensive  and  substantial  smoothing  is  required,  the  use 
of  straight  line  common  tangents  being  a  useful  artifice.  However,  the 
smoothed  velocity  distribution  does  bear  some  resemblance  to  Haselgrove' s 
velocity  distribution.  The  strength  of  the  singularities  associated  with 
this  configuration  gives  a  pointer  that  the  application  of  slender  body  theory 
to  the  effective  aerodynamic  profile  is  likely  to  produce  errors  in  normal 
force  characteristics  and  centre  of  pressure  positions.  Table  1  contains 
these  estimates  and  also  experimental  data  of  Robinson  (17) .  Robinson  used 
a  nose  boundary  layer  trip  and  checked  that  attached  turbulent  flow  occurred 
over  the  whole  body  surface  for  both  bodies  no.  3  and  4.  Slender  body 
theory  is  seen  to  severely  underestimate  the  normal  force  curve  slope  and  to 
predict  a  much  too  far  forward  centre  of  pressure  position.  The  present 
method  predicts  a  normal  force  curve  slope  within  10%  of  the  measured  value, 
which  is  indicative  of  a  reasonably  accurate  estimate  of  boundary  layer 
displacement  thickness  at  the  body  base.  However,  errors  in  estimated 
centre  of  pressure  position  are  large,  particularly  for  body  no.  3  with 
the  smaller  diameter  base.  There  is  thus  no  doubt  that  bodies  no.  3  and  4 
are  aerodynamical ly  non-slender  and  lie  outside  the  range  of  shapes  which 
can  be  treated  by  the  method. 

4.3  Ogive-cylinders 

Ogive-cylinders  can  be  constructed  very  simply  from  two  panels,  which  are 
used  to  represent  the  nose  region.  Figures  5  and  6  show  ogive-cylinders 
with  respective  fineness  ratios  6  and  10  (bodies  no.  5  and  6),  the  S"(x) 
distributions  used,  and  the  comparison  between  the  profiles  achieved  and 
those  desired.  As  indicated  in  Appendix  I  the  "universal"  nose  profile  is 
fatter  than  a  tangent  ogive.  Table  1  contains  estimates  of  the  aerodynamic 
characteristics  based  on  slender  body  theory  and  on  the  present  method. 

The  present  method  has  been  applied  first  of  all  on  the  basis  that  the 
parallel  afterbody  is  a  cylindrical  section,  for  which  the  boundary  layer 
shape  factor  is  1.4,  and  secondly,  by  treating  the  rear  half  of  the  after¬ 
body  as  the  extreme  case  of  a  boattail  of  zero  angle,  over  which,  by  the 
empiricism  in  Section  3,  H  =  3.0. 
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.  Body  3  profile 

.  Body  4  profile 

0 

Profile  produced  by  S"(x)  distribution  A 

L, 

X 

Profile  produced  by  S"(x)  distribution  B 

Parallel  Section 


Boat tail 


S"(x)  distribution  A 
S"(x)  distribution  B 

Haselgrove  (13)  -  Albone's  computer  program 
Slender  body  theory  using  S"(x) 
distribution  A 

Slender  body  theory  using  S"(x) 
distribution  B 

Straight  line  common  tangents 


Extended  boattail-No 
after  body 


(a)  Body  profiles 

(b)  S"(x)  distributions  for  body  no.  3 

(c)  Velocity  distributions 


Figure  4.  Approximation  to  bodies  no.  3  and  4  and  calculated 
velocity  distributions  at  M  =  0,  a  =  0° 
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(a)  Body  profile 

(b)  S"(x)  distribution 


Figure  5.  Approximation  to  body  no.  5 


(a)  Body  profile 

(b)  S"(x)  distribution 

Figure  6.  Approximation  to  body  no.  6 
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The  estimates  for  body  no.  5  are  compared  with  normal  force  data  taken  from 
experiments  by  Ward  (18)  and  an  empirically  determined  centre  of  pressure 
position  given  by  Darling  (19).  Comparative  data  for  body  no.  6  were 
obtained  by  Robinson  (17),  who  used  a  nose  boundary  layer  trip  to  ensure  the 
development  of  a  turbulent  boundary  layer.  For  both  these  shapes  the  boun¬ 
dary  layer  displacement  thickness  is  small  compared  with  the  base  radius  and 
slender  body  theory  gives  close  estimates  of  normal  force  curve  slope  and 
centre  of  pressure  position.  Allowing  for  the  boundary  layer  growth  by 
using  the  present  method  gives  excellent  agreement  with  experiment.  As 
mentioned  in  Section  3,  the  small  changes  resulting  from  treating  the  after¬ 
body  as  a  cylinder  or  as  a  zero-angle  boattail  are  unimportant. 

4.4  Shell  configuration 

The  nose  bluntness  of  body  no.  7  in  figure  7(a)  was  ignored  and  the  body 
represented  by  two  panels  in  the  nose  region  and  one  over  the  boattail 
(figure  7(b).  The  final  entry  in  Table  1  shows  a  comparison  between  results 
of  slender  body  theory,  the  present  method,  and  experimental  data  of  Jermey 
(20) . 


(a)  Body  profile 


(b)  S"(x)  distribution 


Figure  7.  Approximation  to  body  no.  7 
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Jermey's  experiments  were  at  moderately  low  Reynolds  number  and  he  forced 
boundary  layer  transition  on  the  nose  by  a  boundary  layer  trip.  Slender 
body  theory  is  seen  to  underestimate  normal  force  curve  slope  and  to  predict 
the  centre  of  pressure  position  too  far  forward,  whereas  results  using  the 
present  method  compare  closely  with  experiment. 


5.  CONCLUDING  REMARKS 

The  success  or  otherwise  of  the  present  method  appears  to  depend  largely  on 
whether  the  body  being  examined  is  aerodynamical ly  slender.  Slender  body 
theory  defines  what  this  means  mathematically  but  does  not  give  a  quantitative 
explanation.  Criteria  defining  the  limits  of  aerodynamic  slenderness  are  needed 
but  none  are  offered  in  this  report.  Of  the  cases  that  have  been  examined, 
bodies  no.  3  and  4  are  not  aerodynamical ly  slender,  particularly  in  the  vicinity 
of  the  nose,  whereas,  generally  speaking  the  other  bodies  are.  In  any  case 
where  the  body  is  aerodynamical ly  slender  and  the  boundary  layer  displacement 
thickness  is  everywhere  small  relative  to  the  body  radius,  slender  body  theory 
can  be  used  to  give  close  estimates  of  normal  force  and  centre  of  pressure  posi¬ 
tion.  Where  the  boundary  layer  displacement  thickness  becomes  appreciable  in 
comparison  to  the  body  radius,  slender  body  theory  loses  accuracy  and  allowance 
must  be  made  for  boundary  layer  growth.  In  practice,  slender  body  theory  should 
apply  to  configurations  with  favourable,  zero,  or  small  adverse  longitudinal 
pressure  gradients,  such  as  ogive-cylinders  and  ogive-cylinder-flare  combinations, 
provided  they  are  not  excessively  long.  (Esch  (8)  has  shown  that  an  allowance 
for  the  boundary  layer  should  be  made  on  long  bodies) .  When  applied  to  bodies 
with  boattails,  slender  body  theory  seriously  underestimates  the  normal  force 
curve  slope  and  predicts  a  centre  of  pressure  position  which  is  usually  much  too 
far  forward.  In  such  cases  it  is  imperative  that  allowance  be  made  for  the 
boundary  layer  growth.  Predictions  based  on  the  method  in  this  paper  have 
compared  favourably  with  experimental  data  on  an  extensive  range  of  body  shapes 
and  under  conditions  covering  a  large  range  of  boundary  layer  growth  rates. 

The  method  has  proved  versatile  and  reasonably  accurate  for  shapes  which  are 
aerodynamical ly  slender,  giving  the  following  errors: 


Maximum  Error 

Normal  force 

Centre  of  Pressure 

Configuration 

Curve  Slope 

Position 

(%) 

(Calibre) 

Bodies  without 
extensive  boattails 

5 

0.6 

Bodies  with 
extensive  boattails 

3 

2.2 

Increased  accuracy  could  probably  be  achieved  by  allowing  for  boundary  layer 
induced  camber  and  making  a  suitable  compensating  adjustment  to  the  value  of 
the  boundary  layer  shape  factor  in  regions  of  adverse  pressure  gradient,  but 
at  the  cost  of  increased  complexity. 
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NOTATION 

maximum  value  of  S"(x) 

8  * 

boundary  layer  shape  factor  (=  -q  ) 

parameter  related  to  boattail  angle 
body  length 
nose  length 

pitching  moment  about  nose 
free  stream  Mach  number 

UX 

a  Reynolds  no.  (=  —  ) 

body  cross-sectional  area  at  station  x 
cross-sectional  area  at  body  base 


free  stream  speed 

fluid  speed  at  edge  of  boundary  layer 
body  volume 

distance  from  leading  edge  of  equivalent  flat  plate 

minus  normal  force 

constants  in  polynominal  equation 

lengths  defined  in  figure  1 

body  radius 

body  radius  at  end  of  nose  section 

perturbation  velocity 
axial  coordinate 
body  incidence 

=  vl Iff 

O 

a  constant,  whose  value  depends  on  the  magnitude  of  R 

a 


5* 


boundary  layer  displacement  thickness 
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e 
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f 

P 


1,2,3 

B 


boundary  layer  momentum  thickness 

kinematic  viscosity 

a  dummy  variable 

fluid  density 

subscripts 

locations  of  changes  in  body  shape  (figure  8) 
body  base 
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TABLE  1.  COMPARISON  OF  MEASURED  AND  ESTIMATED  NORMAL  FORCE 
CURVE  SLOPE  AND  CENTRE  OF  PRESSURE  POSITION  FOR  MISSILE  BODIES 


Normal  force  slope 

-Cza 

C.P.  position  ahead 
of  nose  (calibres) 

Remarks 

Slender 

Body 

Theory 

Present 

Method 

Experiment 

Slender 

Body 

Theory 

Present 

Method 

Experiment 

Body  No.  1 

Slender 

boattailed 

body 

0.23 

0.77 

0. 79(M  -  0.4) 
o 

0.76  (Me»  0.6) 

37.0 

7.2 

8.81(Mo-  0.4) 

9.39(Mo-  0.6) 

Rl  -  1.9  x  10* 

Boundary  layer  trip 
on  nose 

Experimental  data  ref. IS 

Body  No.  2 

Slender 

boattailed 

body 

0.41 

0.63 

0.65  (Mo-  0.8) 

IS. 8 

8.6 

7.60(Mo-  0.8) 

Rl  •  2.2  x  10* 

Experimental  data  ref. 16 

Body  No.  3 
Non-slender 
boattailed  body 

0.26 

0.S7 

0.63  (M  -  0.6) 
o 

26.1 

10.4 

2.J8(M0.  0.6) 

R,  -  2.4  x  10* 

Boundary  layer  trip 
on  nose 

Experimental  data  ref. 17 

Body  No.  4 

Non-slender 

boattailed 

body  with 

parallel 

afterbody 

O.SS 

0.86 

0.99  (Mo-  0.4) 

0.92  (Mo-  0.6) 

9.b 

S.O 

2.S7(M0-  0.4) 

3.00(Mq-  0.6) 

Rl  -  2.4  X  10* 

Boundary  layor  trip 
on  nose 

Experimental  data  ref. 17 

Body  No.  5 

Ogive- 

cylinder, 

1/d  -  6 

2.00 

2.11 

2.23* 

2.09  (Mq  - 

0  to  0.7) 

-0.92 

-1.01 

-1.20* 

-1.0  to  -1.3 
(M0  *  0.6) 

Rl  •  4.9  X  10* 

•Rear  half  of  afterbody 
treated  as  zero-angle 
boattail 

Normal  force  data  ref. 18 
C.P.  position  ref. 19 

Body  No.  6 
Ogive-cylinder, 

1/d  -  10 

2.00 

2.22 

2.47* 

2.22  (MQ-  0.4) 

2.31  (Mo-  0.7) 

-1.23 

-1.62 

-2.17* 

-2.06 (Mq*  0.S) 

-2.18(M0«  0.7) 

Rl  .  2.1  x  10* 

*Rear  half  of  afterbody 
treated  as  zero-angle 
boattail 

Experimental  djta  ref. 17 

Body  No.  7 
Ogive-cylinder 
with  boattail 

1.47 

1.68 

1.60  (M  •  0.7) 

0.9 

-0.36 

-0.36(Mo*  0.7) 

Rl  >  1.8  x  10* 

Boundary  layer  trip  on 
nose 

Experimental  data  ref. 20 
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APPENDIX  I 

EVALUATION  OF  PERTURBATION  VELOCITIES  ASSOCIATED  WITH  NOSES, 

BOATTAILS  AND  STABILISING  FLARES 

The  three  cases  illustrated  in  figure  1  must  be  examined.  In  each  case  the 
S"(x)  distribution  is  zero  from  x  =  0  until  x  =  f,  and  from  x  =  g  to  the  end  of 
the  body  at  x  =  L.  For  f  <  x  <  g  the  distribution  of  S"(x)  varies  as  shown  in 
figure  1,  in  each  case  attaining  a  maximum  value  of  G. 

I.  Velocity  distribution 

1.1  Constant  S"(x)  distribution 

For  this  case 


f  <  x  <  g  ,  S"(x)  =  G  . 

Inserting  these  values  of  S"(x)  and  S"  (£)  into  equation  (2)  and 
integrating  leads  to  the  following  expressions: 


For 


for 


and  for 


o  <  x  <  f 


u  G  .  f-x 
U0  "  in  g^x  * 


f  <  x  <  g 


iL  =  G  (lr)2 
U  4ir  in  4  (x-f)  (g-x)  ’ 


g  <  x  <  L 


u.  G_  1  x-f 

U  =  4v  ln  x-g  ) 

o  6  •  / 


> 


(1.1) 


1.2  Triangular  S"(x)  distribution  -  increasing  with  x 
For  this  case 


f  <  x  <  g  , 


S"(x)  = 


G  (  x-f 


Again,  inserting  these  values  of  S"(x)  and  S"(£)  into  equation  (2)  and 
integrating  leads  to  the  following  expressions: 


For 


o  ^  x  <  f 


u^ 

U 


4>r  [  _1  +  (  g-f)  ln  g-^  ]  ’ 


(1.2) 
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f  <  X  <  g 


\  g-f 


(g-f)  ln  4(x-(f)(g-x)J  *  (I*2) 


and  for 


g  <  x  <  L 


U0  =  *r  [  U  (g-f)  ln  x-g  ]  • 


1.3  Triangular  S"(x)  distribution  -  decreasing  with  x 


For  this  case 


f  <  x  <  g  S"(x)  = 


;(!3) 


Again,  inserting  these  values  of  S"(x)  and  S"  (£)  into  equation  (2) 
and  integrating  leads  to  the  following  expressions: 

For  o  <  x  <  f 


(S3)  *£] 


f  <  x  <  g 

*  (§3)  *■ 


(1.3) 


and  for 


g  <  x  <  L 


H_  _  G  f  1  .  ( X-g  \  .  x-f  1 

%  '  «L _1  If?)  ln^J  • 


II.  Representation  of  body  components 

The  method  of  employing  the  S"(x)  distributions  will  now  be  described. 
II. 1  Nose  (figure  8(a)) 

For  all  noses  we  assume  that  the  nose  length  is  L  and  that  the 

N 

body  radius  is  rQ  at  x  =  1^. 

II. 1.1  Conical  nose 

For  this  case 

r  =  r  f 

°  hi 


(c)  Parallel  section  -  boat tail  fairing 


Figure  8 

l _ _ _ 


Body  component  parts 


I I. 1.2  Curve  nose 


We  consider  a  curved  nose  which  runs  tangentially  into  the 
next  section  of  the  body  which  is  assumed  to  be  cylindrical, 
and  assume  that,  in  the  context  of  this  paper,  all  such  noses 
can  be  represented  approximately  by  the  equation 


(f.)'  •  i,  •  •(%)••*(  frj 


the  form  of  which  ensures  a  linear  S"(x)  distribution. 
The  boundary  conditions 


v! 


=  i. 


d  r2 
dx 


=  o 


enable  a  and  b  to  be  evaluated  and  lead  to  the  nose  profile 

(1.4) 


(O'  -U-vW] 


This  profile  is  shown  plotted  in  figure  9.  Also  shown  for 
comparison  are  the  nose  shapes  of  several  missile  bodies 
examined  in  Section  4.  It  is  seen  that  the  profile  is  a 
close  representation  of  only  one  body  nose  (that  of  body  no.  1) 
and  it  is  a  fatter  profile  than  a  tangent  ogive.  Furthermore, 
the  profile  is  blunt-nosed,  and  thus  slender  body  theory  cannot 
be  expected  to  hold  near  x  =  o.  Hence,  at  first  glance  the 
chosen  nose  profile  is  not  a  close  approximation  to  many  prac¬ 
tical  nose  shapes.  However,  this  discrepancy  leads  to  very 
small  errors  in  boundary  layer  development  because  the  pressure 
gradient  is  strongly  favourable,  and  the  approximation  is 
adequate . 

Equation  (1.4)  leads  to  an  S"(x)  distribution  given  by 


4 


.  ....  . 


. —  -L  - 


S"  (x) 


H  |  H 
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Figure  9.  Comparison  between  nose  profiles  of  bodies  1  to  7 
and  profile  produced  by  chosen  S"(x)  distribution 


This  can  be  represented  by  the  addition  of  two  panels,  namely 
one  having  constant  S"(x)  distribution  (equation  (1.1))  with 

G  =  2ir  (  —  )  ,  and  one  having  a  triangular  S(x)  distribution  - 

\  N'' 

increasing  with  x  (equation  (1.2)),  with  G  =  -8jt( 

X  hi 

Alternatively,  the  representation  can  be  by  two  different 
panels,  namely  one  having  constant  S"(x)  distribution 
(equation  (1.1)),  with 


and  one  having  a  triangular  S"(x)  distribution  -  decreasing 
with  x  (equation  (1.3)),  with 
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Equations  (1.1),  (1.2)  and  (1.3)  are  evaluated  in  this  case  with 
f  =  0  and  g  =  L^. 

11.2  Boattail  or  flare  (figure  8(b)) 

The  boattail  or  flare  is  defined  such  that  at  a  distance  Li  from  the 
nose  the  radius  is  ri  and  at  distance  L2  the  radius  is  r2 .  Therefore, 
at  any  intermediate  distance  x  the  radius  is  given  by 

r  =  r2  ♦  (n-r2)  (■  .  (1.5). 

The  area  distribution  S(x)  can  readily  be  found  and  thence  S"(x),  given 
by 

s"«  ‘  21  (w)  1  • 

which  is  seen  to  be  the  same  irrespective  of  whether  a  boattail  or 
flare  is  being  considered,  and  is  independent  of  x.  Hence,  the 
contribution  it  produces  to  the  perturbation  velocity  is  given  by 
equation  (1.1),  with  f  =  Li ,  g  -  L*  and 


II. 3  Parallel  portion  -  boattail  fairing  (figure  8(c)) 

For  the  majority  of  missile  body  configurations  a  boattail  is 
preceded  by  a  parallel  section  and  the  junction  between  the  two  regions 
can  be  in  the  form  of  either  a  sharp  kink  or  a  smooth  fairing.  Since 
the  boundary  layer  starts  to  grow  very  rapidly  downstream  of  this 
junction  it  is  desirable  to  obtain  a  fairly  close  approximation  to  the 
velocity  distribution  here,  and  hence  some  attempt  should  be  made  to 
represent  the  fairing  by  an  appropriate  S"(x)  distribution. 

The  fairing  of  length  (Li  -  L3)  is  shown  in  figure  8(c).  It  is 
tangential  to  the  parallel  section  at  (L3,r3)  and  to  the  boattail  at 
(Li,ri).  To  ensure  a  linear  S"(x)  distribution  we  assume  that  the 
radius  squared  can  be  represented  by  the  cubic  equation 

r2  *  a  +  b  (x-L3 )  ♦  c  (x-U)2  +  d  (x-Ij)3  , 

in  which  the  constants  a,  b;  c  and  d  are  found  by  satisfying  the  boun¬ 
dary  conditions 


=  L3  , 

.  d  r2 

r  =  r3  and 

-  0  , 

=  Li  , 

.  d  r2 

r  =  n  and 

=  J  , 
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dr2 

where  J  is  the  value  of  ^  derived  from  equation  (1.5)  when  x  «  Li  , 
namely 


Application  of  the  boundary  conditions  leads  to  the  following  values 
for  the  constants: 


b  =  0, 

‘  -  -  cV,[ 

d  *  anr:>‘  [ 

S"(x)  can  now  be  evaluated  as 

S"(x)  =  2v  c  +  3  d  (x-L3)  j  , 

which  can  be  represented  by  the  addition  of  two  panels,  namely  one 
having  a  constant  S"(x)  distribution  (equation  (1.1))  with 


G  =  2ff  c  , 


and  one  having  a  triangular  S"(x)  distribution  -  increasing  with  x 
(equation  (1.2)),  with 


G  =  6*  d  (L1-L3)  . 

Alternatively,  the  representation  can  be  by  two  different  panels,  one 
having  constant  S"(x)  distribution  (equation  (1.1))  with 

G  =  2*r  (  c  +  3  d  (Li  -L3 )  J  , 

and  the  other  having  a  triangular  S"(x)  distribution  -  decreasing  with 
x  (equation  (1.3)),  with 

G  -  -  6*  d  (L,  -L3 )  . 


Equations  (1 . 1) ,  (1 . 2)  and  (1 . 3)  are  evaluated  with  f  =  Ls  and  g  =  Li . 
The  example  dealing  with  body  no.  1  illustrates  that  the  representation 
of  the  fairing  by  the  cubic  in  r2  is  quite  accurate. 
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The  following  corrections  should  be  made: 

(1)  On  page  22  -  the  middle  equation  of  (1.3)  there  should 
be  a  closing  bracket  |  before  the  comma. 

(2)  On  page  25-7  lines  from  the  bottom,  the  equation 

should  read  G  =  -6n  in  place  of  -8jt 

(3)  On  page  25-3  lines  from  the  bottom,  the  equation 

should  read  G  =  in  place  of 


(4)  On  page  26  -  top  line,  the  equation  should  read 


in  place  of 


sW?5 


These  changes  do  not  affect  any  other  part  of  the 
report . 
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